Effect of shear force on the separation of double stranded DNA 
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Using the Langevin Dynamics simulation, we have studied the effects of the shear force on the 
rupture of short double stranded DNA at different temperatures. We show that the rupture force 
increases linearly with the chain length and approaches to the asymptotic value in accordance 
with the experiment. The qualitative nature of these curves almost remains same for different 
temperatures but with a shift in the force. We observe three different regimes in the extension of 
covalent bonds (back bone) under the shear force. 
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Inter- and intra- molecular forces are key to the sta- 
bility of DNA and biological processes e.g. transcrip- 
tion, replication, slippage etc. [HQ. Up to now, under- 
standing of these forces was possible through the indirect 
physical and thermodynamical measurements like crys- 
tallography, light scattering, nuclear magnetic resonance 
spectroscopy etc. Q. Single molecule force spectroscopy 
(SMFS) experiments have directly measured these forces 
and provided unexpected insights into the strength of the 
forces driving these biological processes as well as deter- 
mined various interactions responsible for the mechani- 
cal stability of DNA structures 043- With the increas- 
ing number of experiments and insights gathered so far, 
it has become clear that the measurement of molecular 
interactions not only depends on the magnitude of the 
applied force, but also depends how and where the force 
was applied |7H15|. 

A major concern is now to understand whether all 
these interactions contribute at the same moment or they 
have different life times. In order to understand this, 
a force has been applied perpendicular to the helix di- 
rection (DNA unzipping) and along the helix direction 
(rupture and slippage) as shown in Fig. 1 043 ■ In- 
case of unzipping of double stranded DNA (dsDNA), the 
critical force is found to be independent of the length 
of DNA and the loading rate @,|^). This may be un- 
derstood theoretically that at the centre point of fork 
(Fig. lb), the applied force only breaks a base pair at 
a time and hence it remains independent of the loading 
rate and length. However, when a force (up to 65 pN) is 
applied along the helix direction (shear force) , the length 
of the dsDNA increases and the force-extension (/ — x) 
curve can be described by the worm like chain (WLC) 
model [I?}. In the high force regime (> 65 pN), the 
dsDNA can be overstretched about 1.7 times of the B- 
form contour length and a phase transition occurs from 
the B-form to a stretched or S-form [l8l - [20j . Recently, 
van Mameren et al. have studied DNA stretching with 
or without DNA binding ligands and demonstrated that 
overstretching comprises a gradual conversion from ds- 
DNA to ssDNA and it should be interpreted in terms of 
force induced DNA melting [2~l| . 

For a short dsDNA, if the applied shear force increases, 
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FIG. 1: Schematic representation of dsDNA: (a) dsDNA in 
zipped form; (b) Unzipping of dsDNA by the force (/) applied 
at one end (5' — 3'); (c and d) Shear force along the chain 
applied at the opposite ends (3' — 3' or 5' — 5') of the dsDNA. 



the dsDNA separates into two single strands at some crit- 
ical force. This phenomenon has been identified as rup- 
ture li"oj . The unbinding force strongly depends on 
the pulling end and is much larger than the unzipping 
force [1 Hfl, Hill- Neher and Gerland studied the dy- 
namics of dissociation of the two strands and found the 
expression for the critical force [23[ . Expressing the bond 
energy and the base pairing energy in the form of har- 
monic oscillators in the ladder model of dsDNA of length 
L, de Gennes [24j proposed the maximum force required 
for the rupture 
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Here, fx is the force required to separate a single 
base pair, which is same for the homo-sequence and 
X^ 1 = \JQ/2R is the de Gennes characteristic length 
over which differential force is distributed. Here, Q and 
R are the spring constants, characteristic of stretching of 
backbone and hydrogen bonds, respectively. 

Recently, Danilowicz et al. [l3[ systematically stud- 
ied the DNA rupture by varying the length of dsDNA. 
The critical shear force is found to increase linearly up 
to a certain length and approaches the asymptotic value 
(s=s 62 pN), which is in good agreement with the de 
Gennes prediction. It was argued that the covalent bonds 
(backbone) and the hydrogen bonds involved in the base 



pairing will be stretched under the applied force. The 
differential force will approach to zero at the length x 1 ^ 
if one moves in from the either side. However, no ex- 
perimental effort has been made to study the effect of 
shearing force on the stretching of covalent bonds and 
hydrogen bonds in side the characteristic lengthy -1 . 
Moreover, in the description of de Gennes model [24| or 
subsequently improved model by Chakrabarti and Nel- 
son [23, effect of thermal fluctuation has been ignored, 
whereas all the rupture experiments were generally per- 
formed at finite temperature. The aim of this manuscript 
is to study the effect of temperature (T) on the rupture 
and consequences of differential force on the distribution 
of extension in bond lengths and hydrogen bonds near 
the rupture. 

We use Langevin Dynamics (LD)simulation to inves- 
tigate mechanical and physical properties related to the 
rupture of DNA [H, H3, [29|, Since, rupture time is 
of the order of milliseconds to seconds, an atomistic sim- 
ulation of longer chain in the solvent is computationally 
difficul t f3lL [331 . We have used a coarse-grained model 
0, [H, l30l. l33j of the flexible polymer chain to model a 
DNA, which allows us to study a larger system size and 
events of a longer time scale. A chain in the model con- 
sists of bead units connected by effective bonds charac- 
terized by the stiff springs. Each effective bond repre- 
sents several chemical bonds (e.g. sugar phosphate etc. 
) along the chain backbone. The energy of the model 
system is given by 
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where N is the number of beads in each strand. iD 
represents the position of i th bead on I th strand. In 
present case, I = 1(2) corresponds to first (compli- 
mentary) strand of dsDNA. The distance between intra 
strand beads, ufj , is defined as \i£p —u^p \ . The harmonic 
(first) term with spring constant k (=100) couples the ad- 
jacent beads along the two strands. Second term takes 
care of excluded volume effect i.e. two beads can not 
occupy the same space [34]. The third term, described 
by Lennard- Jones (LJ) potential, takes care of the mu- 
tual interaction between two strands. The first term of 
LJ potential (same as second term of Eq.2) will not al- 
low the overlap of two strands. Here, we set C = 1 and 
A = 1. The second term of LJ potential corresponds to 
the base pairing between two strands. The base pairing 
interaction is restricted to the native contacts (&y = I) 
only i.e. i th base of I s * strand forms pair with the i th 
base of 2 nd strand only as shown in Fig. la. This is 
similar to the Go model [28]. The parameter do(= 1.12) 
corresponds to the equilibrium distance in the harmonic 
potential, which is close to the equilibrium position of 
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FIG. 2: (a) Force vs extension curves for different chain 
lengths Arrows indicate the maximum force, where the num- 
ber of contacts approaches to zero. For the sake of compari- 
son, we have normalized the extension by its contour length 
(at f = 0). (b) Figure shows the variation of rupture force 
with length. The solid line corresponds to a fit of Eq. 1. 
Solid circles represent the value obtained through the simula- 
tion. 



the average LJ potential. In Eq. 2, we use dimensionless 
distances and energy parameters. The major advantage 
of this model is that the ground state energy of the sys- 
tem is known [28| . The equation of motion is obtained 
from the following Langevin equation [26|, [2?], H3] 
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where m and £ are the mass of a bead and the fric- 
tion coefficient, respectively. Here, F c is defined as 
— 4^ and the random force F is a white noise [27|, he., 
< r(i)r(i') >= 2(TS(t-t'). The choice of this dynamics 
keeps T constant throughout the simulation for a given 
/. The equation of motion is integrated by using the 
Qth orc [ er predictor-corrector algorithm with time step 
St=0.025 [27[. The results are averaged over many trajec- 
tories. The equilibration has been checked by monitoring 
the stability of data against at least ten times longer run. 
We have used 2 x 10 9 time steps out of which first 5 x 10 8 
steps have not been taken in the averaging. 

In the constant force ensemble, we add an energy —f.x 
to the total energy of the system given by Eq. 2. We cal- 
culate the reaction coordinate x (extension) for different 
values of /. The f — x curves (Fig. 2a) show the entropic 
response at low forces and remai n q ualitatively similar to 
the one seen in experiments [f| [lfl EH, E3 • We identify 
the rupture force as a maximum force, where the number 
of intact base pairs suddenly goes to zero. In Fig. 2b, we 
show the rupture force as a function of the chain length 
at low temperature. It is evident from this plot that the 
rupture force approaches to an asymptotic value for the 
chain length greater than 20, which is in accordance with 
the experiment [13j . 

We expand the LJ potential given in Eq.2 around its 
equilibrium value. The coefficient of second (harmonic) 
term of its expansion corresponds to the elastic constant 
of the base-pairing. The de Gennes characteristic length 
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FIG. 3: (a) Variation of extension in hydrogen bond length 

(A h ) along the chain at three different forces, (b) Figure FIG. 4: (a) Variation of the rupture force with chain length 
shows the variation in extension of covalent bond length (A c ) at different T. (b) Figure shows the DNA rupture force- 
along the chain. Here, open and solid symbols correspond to temperature diagram for chain length 32. 
one strand and its complementary strand, respectively. 



[24j for the present model is estimated to be ~ 10. Sub- 
stituting the value of /i(= 1) and the above mentioned 
value of x _1 m Eq. 1, we obtained the value of f c for 
a given length of dsDNA, which is shown by solid line 
in Fig 2b. One can notice a nice agreement between the 
simulation and the value predicted by Eq. 1 [2J |. 

One of the important findings of the present simula- 
tion is the distribution of stretching of hydrogen bonds 
(A^) and extension in the covalent bonds (A c ) for a wide 
range of force below the rupture, which are experimen- 
tally difficult to obtain. In Fig. 3a, we depict the vari- 
ation of Aft with base position for the chain of length 
40. From this plot, one can observe that the hydrogen 
bonds at extreme ends (up to « 10 bases) get stretched, 
whereas bases in the middle (above the de Gennes length 
~ 10 ~ 30) remain same indicating that the differential 
shear force approaches to zero in this region. In Fig. 3b, 
we show the variation of A c with the base position. All 
curves have three distinctively different regions. One can 
observe that bonds near the pulling end (say 5 '-end) get 
stretched more and decreases gradually. However, when 
one approaches the other end (i.e. 3'-end), there is a 
change in the slope and the extension is quite less com- 
pare to the middle one. It should be noted that 3 '-end 
is near to 5 '-end of the other chain, where a similar force 
is also applied. Since, dsDNA is in the zipped state, the 
applied force at 5 '-end of one strand also pulls the other 
strand along the opposite direction, which causes a rela- 
tively slower increase. 

In model studies either temperature is set to zero [24| 
or thermal fluctuation is ignored [25| as a result, the rup- 
ture force defined in Eq. f is independent of temperature. 
However, all rupture experiments usually performed at 
room temperature [f(| E3|, therefore, it is desirable to un- 
derstand role of entropy, which may play a significant role 
at higher T. The thermodynamics of force induced melt- 
ing can be obtained from the following relation [35|, HH : 

-fx = AH- TAS, (4) 

where H is the enthalpy and S is entropy of the ruptured 
chains. Setting x equal to unity and replacing the value 



of AH by the value of the rupture force at T = 0, Eq. 4 
can be written as 

-f = f c - TAS, (5) 

where f c is given by Eq. I. In the thermodynamic limit, 
value of S may be analytically estimated [34| or can be 
obtained from the experiments [37]. However, for the 
finite chain length, one has to resort on numerical tech- 
niques to get the value of S. It is possible to study the 
effect of temperature on the rupture force in the present 
setup. In Fig. 4a, we show the dependence of the rupture 
force on the length at different temperatures. The quali- 
tative nature of the curves obtained at different T remain 
similar to de Gennes plot (Fig. 2b), with a shift showing 
that the rupture force decreases with T. From Fig. 4a, 
one can notice that for chain length 32, the rupture force 
has approached its asymptotic value for all T. In Fig. 4b, 
we depict the force-temperature diagram for the DNA 
rupture for the chain length 32. A linear dependence on 
temperature can be noticed, which is in accordance with 
Eq. 5. Usually rupture experiments are performed well 
below the melting temperature, which correspond that 
the DNA is in the zipped state. Therefore, in Fig. 4b, 
maximum temperature has been set slightly below the 
melting temperature i.e. T m = 0.23 at / = 0. 

In this paper, we have studied the effect of shear force 
and temperature on the rupture of dsDNA. We show that 
the shear force increases linearly with the length of the 
DNA and then approaches to the asymptotic value. In 
the lattice model, the bond length is constant (stiff), 
therefore, covalent bonds and hydrogen bonds will not 
be stretched. As a result, the rupture force increases 
linearly with length (36| in the lattice model or in the 
models, where bond length is considered as a constant. 
Our simulations confirm that, one will not gain strength 
by increasing a larger pair sequence as predicted by de 
Gennes [241. Interestingly, recent experiment also sup- 
ports it [13j . The distribution of extension in hydrogen 
bonds and covalent bonds for different forces are shown 
in Fig. 3, where Fig. 3a clearly shows that the increase 
in the extension of the hydrogen bonds is limited up to 
the de Gennes length, which is consistent with the strain 
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profile obtained by Chakrabarti and Nelson (25{. Above 
this length, differential shear force approaches to zero, 
as a result there is no extension in the hydrogen bonds. 
We also find that the qualitative nature of dependence of 
rupture force on length remains same for different tem- 
peratures with a shift. The rupture force decreases with 
temperature (Fig. 4b) as predicted by Eq. 5. 

Our study shows that the de Gennes length remains 
independent of the applied force. The most surprising 
finding of the present simulation is revealed from Fig. 3b, 
which shows variation of extension in the covalent bonds 
along the chain for three different forces (/ = 10, 14, 18). 
Although, in all these cases, the differential force ap- 
proaches to zero and there is no relative increase in the 
bond length above the de Gennes length along the chain 
(Fig. 3b). However, unlike the hydrogen bonds, there is a 
net increase in the extension in the covalent bonds, which 
depends linearly on the applied force. Nuclear magnetic 
resonance experiment [38| or the atomistic simulation [32[ 
should be able to observe this. 

It may be noted that the present simulation is carried 
out in the reduced unit. It is possible to extract a rough 
estimate of the rupture force in the real unit. The free 
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